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NEARLY RADIAL NEUMANN EIGENFUNCTIONS ON SYMMETRIC DOMAINS 


NILIMA NIGAM, BARTLOMIEJ SIUDEJA AND BENJAMIN YOUNG 


Abstract. We study the existence of Neumann eigenfunctions which do not change sign on the 
boundary of some special domains. We show that eigenfunctions which are strictly positive on the 
boundary exist on regular polygons with at least 5 sides, while on equilateral triangles and cubes it 
is not even possible to hnd an eigenfunction which is nonnegative on the boundary. 

We use analytic methods combined with symmetry arguments to prove the result for polygons 
with six or more sides. The case for the regular pentagon is harder. We develop a validated nu¬ 
merical method to prove this case, which involves iteratively bounding eigenvalues for a sequence 
of subdomains of the triangle. We use a learning algorithm to hnd and optimize this sequence of 
subdomains, making it straightforward to check our computations with standard software. 


1. Introduction 

We study the existence of an eigenfunction of the Neumann Laplacian which is positive (or 
nonnegative) on the boundary of highly symmetric domains. Recently, Hoffmann-Ostenhof [18] 
proved that on rectangles, any Neumann eigenfunction that is positive on the boundary must be 
constant. In this paper we prove similar results for regular polygons and higher dimensional boxes. 

Schiffer’s conjecture (see [38]) states that if a Neumann eigenfunction is constant on the bound¬ 
ary of a domain, then either the eigenfunction is constant in the domain, or the domain must be a 
disk. This conjecture is still open, although many partial results are known (see e.g. [8, 9, 14]). 
We relax the boundary restriction (positive instead of constant) and ask if the modified conjecture 
holds for a special class of domains. 

Our problem has a rather interesting physical interpretation in terms of a sloshing liquid in a 
cup with a uniform, highly symmetric cross-section (see [21] for a relation between sloshing and 
Neumann eigenvalue problem). It is nearly obvious that one can disturb a fluid in a round cup so 
that the created wave is radial. In particular the fluid level rises and lowers simultaneously along 
the whole cup wall. Hoffmann-Ostenhof [18] proved that it is possible to create a wave in a square 
cup so that there are a few stationary points along the wall, but it is impossible to make all points 
move in unison. We prove that no such wave can be created in a triangular cup, even if stationary 
points are allowed. At the same time, it is possible to create a wave with unison movement along 
the boundary of regular polygons with at least 5 sides. In summary, we have: 

Theorem 1.1. Any Neumann eigenfunction that is nonnegative on the boundary of an equilateral 
triangle is constant inside. 

Theorem 1.2 (Hoffmann-Ostenhof [18]). Any Neumann eigenfunction that is positive on the bound¬ 
ary of a rectangle is constant inside. 

Theorem 1.3. There exists a Neumann eigenfunction on a regular polygon with n > 5 sides that 
is positive on the boundary and not constant. 
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Remark. Squares are in some sense a critical case for regular polygons. An eigenfunction that is 
positive on the boundary does not exist, yet 

(p{x, y) = — cos 71X — cos Tiy 

is an eigenfunction of the square [—1,1]^. It is positive on the boundary, except at the midpoints of 
all sides (where it equals 0). 

We also study higher dimensional boxes. Surprisingly, eubes no longer have nonnegative eigen- 
funetions. 

Theorem 1.4. Any Neumann eigenfunction that is nonnegative on the boundary of a cube (more 
generally a box) in dimension d > 2 is constant inside. 

Careful Finite Element computations suggest the following conjecture. 

Conjecture 1.5. Eigenfunctions which are nonnegative on the boundary of a tetrahedron and 
octahedron do not exist. However, eigenfunctions which are positive on the boundary exist on 
dodecahedron and icosahedron. 

The paper uses a variety of methods to handle the different cases. In particular, we use com¬ 
binatorial and number theoretic results on cubes (Section 4) and equilateral triangles (Section 3). 
We dissect regular polygons with n > 5 into congruent right triangles and study their Neumann 
eigenfunctions. For n > 6 we can use existing results on the shape of the second Neumann eigen¬ 
function to draw the necessary conclusions (Section 5). 

The proof for the regular pentagon is rather unusual and the most interesting part of the paper. 
We prove that the nodal line for the second Neumann eigenfunction of a right triangle must connect 
two longest sides. This seemingly simple fact is extremely hard to prove. Similar results for 
obtuse triangles have been obtained by Atar and Burdzy [5] using very sophisticated probabilistic 
techniques. 

Our proof uses a generalization of the recently-described computable lower bounds for eigen¬ 
value approximations [11] in the setting of validated numerics, see Section 8.2. We find lower 
bounds for mixed eigenvalues of 20 complicated polygonal domains and use these in a iterative 
procedure to restrict the shape of the eigenfunction. The proof itself is human-readable and fully 
analytic (Section 6), except for the matrix eigenvalue computations on 20 large matrices from 
Section 8.4 and Section 8.5. We employ three different methods (intentional redundancy to ensure 
correctness, using different software packages) to find lower bounds for the smallest eigenval¬ 
ues of these matrices: LDL' decomposition with interval arithmetic, LU decomposition using 
exact rational representations, and Hessenberg decomposition with interval arithmetic and Sturm 
sequences. 

The 20 domains mentioned above were found using a non-validated learning algorithm quickly 
examining thousands of cases and producing a sequence of domains for the main validated numer¬ 
ical algorithm. Both procedures use a novel approach which iteratively constrains the location of 
the nodal line for the second Neumann eigenvalue in the original triangular domain. Section 7. Our 
learning algorithm is not unique in that one could design a different strategy resulting in a different 
set of the domains, leading to a different validated proof. 
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2. Definitions and auxiliary results 


The Neumann eigenvalue problem ean be approaehed elassieally, by solving the partial differ¬ 
ential equation 

in ^ i 

dyUn = 0 on dD. 


However, it is often more useful to work with the variational weak formulation 




inf sup 

Sc wi-2(d) ues 

dim S = n 




( 1 ) 


where H^{D) is the Sobolev space of all functions u G L‘^{D) such that Vu G {L‘^{D))‘^. The right 
side of (1) is commonly called the Rayleigh-Ritz quotient. In this context the minimizers of the 
Rayleigh quotient are the eigenfunctions. For Lipschitz domains (and even more general domains 
for which appropriate Sobolev embeddings exist) the two approaches lead to the same eigenvalues 
and the same eigenfunctions (via elliptic regularity considerations). For a broad overview on this 
topic see [7] and [10]. 

Note that the variational characterization lacks any obvious boundary conditions. This is a con¬ 
sequence of the Neumann (also called natural) boundary condition being automatically enforced 
by the Sobolev spaces. In contrast, to enforce Dirichlet boundary condition one seeks minimizers 
of the Rayleigh-Ritz quotient over a subset of H^{D) consisting of functions with zero trace on 
the appropriate part of the boundary of the domain. 

In general it is true that 

0 = /Xi < /r2 < /^s < • • • < oo. 

However in some special cases one can show that /r 2 is simple. In particular. 


Lemma 2.1 ([33, Theorem 1]). Fornon-equilateral triangles ^2 simple. 


We will also work with a mixed Dirichlet-Neumann eigenvalue problem: 

\yiIIji in 77, 

= 0 on 77 C dD, 
duUn = 0 on dD \ B. 

Note that eigenfunctions satisfy Dirichlet boundary conditions on 77, and the appropriate varia¬ 
tional formulation must include the same restriction. 

Xn= ml sup— 5 —, 

S(1H],[D)u€S 

dim S=n 

where H^{D) is a subspace of H^{D) consisting of all functions satisfying u = 0 on 77. If 
meas{B) > 0, we have 

0 < Ai < A2 ^ ^ Xn —^ 00. 

In what follows we need a few geometric results from [33]. 


Lemma 2.2 ( [33, Lemma 4] ). Suppose D is a domain with a line of symmetry. Then there cannot 
be two orthogonal antisymmetric eigenfunctions in the span of the eigenspaces 0 //X 2 and (note 
that P 2 might equal ps). 
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Lemma 2.3 ( Special case of [33, Lemma 5]). The nodal line for the second Neumann eigenfunc¬ 
tion on a triangle must start on one side and end on another side or vertex connecting the other 
two sides. 

Let us introduce the following naming convention for isosceles triangles: 

Definition. A triangle is superequilateral (subequilateral) if it is isosceles with aperture angle 
larger (smaller) than tt/S. 

Extensive numerical studies suggest that Lemma 2.3 can be strengthened to: 

Conjecture 2.4. The nodal line for the second Neumann eigenfunction ends in a vertex only for 
superequilateral triangles. For all other non-equilateral triangles the nodal line connects two 
longest sides. 

We prove this conjecture for some right triangles: namely, the ones whose smallest angle is 
approximately vr/S. Note that an even stronger conjecture was posed by Atar and Burdzy [5, 
Conjecture 3.2] for obtuse triangles, where the nodal line would be confined in the right triangle 
bounded by two long sides and the altitude perpendicular to the longest side. 

3. Equilateral triangle: proof of Theorem 1.1 

The Neumann spectrum for an equilateral triangle can be split into symmetric modes pm,n and 
antisymmetric modes forming eigenspaces of eigenvalues Xm,n, n > m > 0. Note that Xm,n 
might be equal for different pairs of integers m, n. This means that there exist eigenfunctions of 
an equilateral triangle that combine many different modes. Eor more details see McCartin [27]. 

Theorem 8.1 from [27] states that the symmetric modes never vanish, while the antisymmetric 
modes degenerate when m = n. Eurthermore, Theorem 8.2 from [27] gives that both types of 
modes are rotationally symmetric if and only if m = n (mod 3). Einally, any symmetric mode 
that is not rotationally symmetric can be written as a sum of two rotated (by 120 and 240 degrees) 
antisymmetric modes. More precisely, denoting the rotations of '<pm,n by 120 and '0m,n,240 we 

have Prn,n 0m,n, 120 "h 0m,n, 240- 

Suppose / is an eigenfunction for some A for an equilateral triangle with horizontal side Si. 
Then / is a linear combination of the symmetric and antisymmetric modes (with respect to the 
altitude Oi perpendicular to Si). If (prn,n is symmetric but not rotationally symmetric, we rewrite it 
using two antisymmetric modes. Therefore 

f= E ^ ^ 120 '0m,n,24O ) + E 

m,k>0 m^n (mod 3) m^n 

Suppose / is nonnegative on the boundary; let G be the group of isometries of the equilateral 
triangle. Then f o U, U G G is also nonnegative on the boundary. Eurthermore, the orbit of G on 
any antisymmetric mode 0 has size 6 (all possible rotations and reflections are different) or length 
2 (if 0 is rotationally symmetric). At the same time, the sum of the reflections of 0 along the line 
of antisymmetry, and in particular at the midpoints of the boundary edges, cancel out. Therefore 
the function 

^ ^ ^ ° G 6 ^ ^ (^m,m+3kTm,m+3k- (2) 

UgG m,k>0 

is also nonnegative on the boundary. It follows that Theorem 1.1 needs to be proved only for 
eigenfunctions of the form F. 
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Consider the equilateral triangle with vertiees (0, 0), (1, 0) and (1/2, \/3/2). We have the fol¬ 
lowing symmetric modes (see e.g. [27]) 


( 3 ) 


y) = COS ( -7r(2x - l)(m - n)j cos + n)y ] + 

+ (—1)”^ cos ( ^vr(2a: — 1)(m -|- 2n) j cos ( -^nmy ) -|- 


-h (—1)”' COS ( -vr(2a: — l)(2m -f n) j cos ( —^nmy ) , 


v/3 

2 

7 !' 


with corresponding eigenvalues 


Am.r). 




(m^ -|- mn + n^). 


Note that an eigenvalue A might have a high multiplicity, since many different pairs of integers 
(m, n) might produce the same value A. Therefore more than one pair of integers might belong to 
a particular eigenvalue A. 

We find that the rotationally symmetric modes satisfy 


0) = cos{27ikx) + cos(27r(m + 2k)x) -f cos(27r(m + k)x), on [0,1]. (4) 


Note that we only need to consider the values of ^rn,m+ 3 k on one side, so we assumed y = 0. This 
last quantity integrates to 0 on (0,1), the side of the triangle, unless k = 0. This means that a pair 
(M, M) must belong to A if the eigenfunction associated to A has nonnegative boundary values. 
This implies that all other pairs (m, n) for the same A satisfy 3M^ = + mn + n^. 


3.1. Eigenfunctions positive on the boundary. Note that if both m and n are even, then m? + 
mn + is even, otherwise it is odd. Therefore an even M implies even m, n. By induction, any 
pair (m, n) belonging to eigenvalue A must have both m and n divisible by 2® and at least one not 
divisible by 2®+^, whenever 4® divides M^, but 2 • 4® does not. Hence there exists s such that for 
any pair (m, m + 3A:) that belongs to A we must have m = 2^mi and k = 2^ki, where at most one 
of the ki and mi is even. Therefore 

</3m,m+3fc(2“''“\ 0) = cos(7rA;i) + cos(7r(mi + 2A;i)) + cos(7r(mi + ki)) = 

= (- 1)^1 + (- 1 )”^! + = - 1 , 

^m,m+3fc(0) 0) 3. 

We have proved that all (pm,m+ 3 k that belongs to A must equal —1 at the same point on the boundary 
(and 3 at vertices). Hence any linear combination of such eigenfunctions with ^ am,m+ 3 k 7^ 0 
in F must change sign. Also, unconditionally the eigenfunction cannot be strictly positive. 

This argument is very similar to the one used by Hoffmann-Ostenhof [18] on squares to prove 
nonexistence of eigenfunctions positive on the boundary. It is however impossible rule out the 
existence of eigenfunctions nonnegative on the boundary using this method. We might be able to 
prove that the linear combination must equal 0 at many points, but not that it changes sign. 


3.2. Eigenfunctions nonnegative on the boundary. Here we develop an improved method based 
on the fact that Neumann eigenfunction cannot vanish on an open subset of the boundary (it would 
satisfy Dirichlet condition). Therefore, if we find an open set on which eigenfunction integrates to 
0 , it must change sign on that set. 
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We ony need to work with A = (m = M, fc = 0 is admissible in (4), see eomment 

below that formula). It is also possible that other pairs (m, k) with k > 0 give the same A. In that 
ease we have 

3M^ = 3m^ + 3mk + k"^. 


It is easy to eheek that 

M<m + k<m + 2k< 2M, 0 < k < M. 


( 5 ) 


Consider points 


2i + l 

Xi = - 

2M 


0<i< M, 


and integrals over symmetrie intervals around these points 


M—1 pxi—a 2 M—1 

/ cos(aa;) dx = — sin(Q!a) cos(aa;j) 

i=0 Jxi+a ^ 


— sin(aa) cos(a;/2) sin(a/2) csc(a/2M) 
a 


by [16, Seetion 1.341, Formula 3], as long as sm{a/2M) ^ 0. 

Examining formula (4) we find a = 2'Kk, 27r(m+A:) and 27r(m+2fc). In eaeh ease sin(a/2M) ^ 
0 due to (5). Furthermore sin(Q;/2) = 0, henee the eigenfunetions (^m,m+ 3 fc with fc > 0 integrate 
to 0 over the union of (xj — a^Xi + a). We only need to show the same property for for the 
eigenfunetion with fc = 0: 


0) = 1 + 2 cos{27iMx). 


( 6 ) 


We have 


M—l pxi—o 


M-l 


j=0 


j=0 


I ipM,M{x,0) dx = y2a -\—^ sin(27rMa) cos((2i + l)7r) 


= 2aM -sin(27rMa). 

TT 


Let zq = 271 Mao and find a positive solution of zq = 2sin(2;o)- We get zq < ti and qq < 1/2M. 
Henee intervals (x* — oq, x* + oq) fit inside (0,1), the side of the equilateral triangle. At the same 
time, any linear eombination of eigenfunetions from (4) integrates to 0 over the union of these 
intervals. Henee it must ehange sign, as it eannot satisfy both Diriehelt and Neumann eonditon on 
any interval (be identieally 0 on any interval). 


4. Cubes: proof of Theorem 1.4. 

4.1. Fully symmetric eigenfunctions. Consider the eube C = [—1,1]"^. Suppose it has a Neu¬ 
mann eigenfunetion that is positive (nonnegative) on the boundary. We ean symmetrize this eigen¬ 
funetion by applying all isometries of the eube and summing the resulting eigenfunetions (as in 
the equilateral triangle ease). We will get a new eigenfunetion that is positive (nonnegative) on the 
boundary, symmetrie with respeet to Xj = 0 for any i and invariant under arbitrary permutation 
of variables x^. We only need to prove that this fully symmetrie eigenfunetion eannot be positive 
(nonnegative) on the boundary. 

Any symmetrie eigenfunetion ean be written as a sum of simple eigenfunetions of the form 

n 

(_l)Emi JJ" cos(mi7rxj). 
i=l 
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The factor (—ensures positivity in all vertices (Xj = ±1). Invariance under permutations of 
variables gives 




Vm? = A i=l 

^ i 

M={m\ <---<mn} 


(7) 


where ct„ denotes any permutation of {1,..., n} and qm are arbitrary coefficients depending on the 
nondecreasing sequence of nonnegative integers rrii. We require that ^ = A to ensure all terms 

belong to the same eigenvalue. Formula (7) gives the most general form of the eigenfunction that is 
invariant under the group of the isometries of the cube C. We need to show that this eigenfunction 
is negative somewhere on the boundary of the cube, regardless of the choice of A. Due to symmetry 
we only need to consider one face. 

Note also, that ^px is also a linear combination of eigenfunctions of the lower dimensional Lapla- 
cian on a face. Indeed, fixing xi = 1 gives a sum of products of cosines, hence again a symmetric 
function. However, due the presence of the permutations an, we drop different in different 
terms, and hence we get a sum of eigenfunctions for various eigenvalues. Every non-constant 
Neumann eigenfunction is orthogonal to the constant eigenfunction. Hence ipx integrates to 0 on 
each face, unless an eigenfunction which is constant on the face is a part of ipx, cf. the discussion 
below (4) pertaining to equilateral triangles. 

Therefore the sequence mi = ■ ■ ■ = mn-i = 0, = m = \/A gives one of the terms in ipx. 

Consequently, \ = m? for some integer m. Otherwise ipx integrates to 0 over any face, hence it 
must change sign on each face. 


4.2. Positive eigenfunctions. We begin with a special case to illustrate the approach. Suppose 
X = m? with odd m. Recall that ipx is a sum over all sequences mi < m 2 < ■ ■ ■ < such that 


2 I I 2 \ 2 

m^^ + • • • + = X = m . 


Hence at least one is odd. Consider a discrete set of points: 

= {(a:i > X2 > ■ ■ ■ > Xn) ■■ Xi e {Q, 1}}. 


These points correspond to the center of the cube (0,..., 0), the center of the face (1,0,..., 0), 
the centers of all lower dimensional faces, finally a vertex (1,..., 1). Let Xq be the set 


Xq = {{1 = Xi > X2 > ■ ■ ■ > Xn = 0) : Xi e {0, 1}}. 
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Note that all points in Xq are on one faee of the eube. For any x E X put k = (*^he 

eodimension of the faee for whieh x is a eenter). Observe that 


x^Xq 


n—\ 


1 

Ed 

1 

1 

[n — k)\ 

k 

■ 5^ n cos("ia„(i)vr 

k=l 



an i=l 

n—1 

= E 

M={mi <"-<mn} 





1 

1 

{n — k)\ 

■ y^( X)55i-1 

k=l 

Em|=> 


<Tn 


M={mi <-"<mn} 

71—1 

55 

k=l 



= 

aM(-l) 

Em2=A 

1 

(n — k) 

- ^ ^( ^an(i) 

O'n 




Note that in the innermost sum eaeh term appears exaetly {n — k)\ times, sinee we are using only 
the first k values of eaeh an- Henee we are adding (exaetly onee) all produets of (—1)'"% exeept 
for the full and empty produet, so we may rewrite this as 

^ % 

But at least one is odd, henee the produet in the braeket is 0. Furthermore, the sum of m* is 
also odd, henee the whole bracket is 0. Thus 


JJ(1 + (-1)™*) - (-1)55”^* _ 1 

i=\ 


5 : 

x&Xq 


{n 


1 


^x{.x) = 0 . 


Therefore either ipx is 0 at the centers of faces of arbitrary dimension, or ifx must change sign. 
To prove the eigenfunction must change sign we will use a different method, similar to the one 
for equilateral triangles (Section 3.2). For the moment, we can show that eigenfunction cannot be 
positive on the boundary for a few low-dimensional cases with an argument about the parity of the 

rrii. 


Proposition 4.1. In dimensions 2, 3 and 4, any positive Neumann eigenfunction on a cube must be 
constant. 


Remark. Dimension 2 was proved by Hoffman-0 stenhof\l%]. 

Proof. We only need to consider even m. For 0 < h < 1 define 

Xh = {{I = Xi> X2 > ■ ■ ■ > Xn = h) -. Xi e {h, 1}}. 

As above we get 

n 

]^(cos(mi7rh) + (—1)”**) — 1 — J^cos(mj7rh) , 

. 2=1 i 

M={rni<---<rriri} 
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We now consider each dimension separately: 

• Dimension 2: 

The sum of the squares of two odd numbers is congruent to 2 modulo 4, hence it is not 
a square. Therefore, if m is even, then both mt are even. Furthermore, by induction 2^ 
divides both mj, but 2®+^ divides exactly one of them. Take h = 1/2®. Then cos(mj7r/2®) 
have both signs. But both irii are even, hence the first product in the bracket is 0, and the 
second product equals —1. Hence the whole bracket equals 0. 

• Dimension 3: 

The sum of three squares is again a square only if all numbers are even. Indeed, with 
two odd numbers, the sum that is congruent to 2 modulo 4. By induction, there exists s 
such that 2® divides all m*, while 2®+^ divides none or two. In either case, the first product 
equals 0, and the second equals —1. Hence the bracket is again 0. 

• Dimension 4: 

The sum of k odd squares is congruent to k modulo 8. Hence only 1 or 4 odd squares 
can give a square. Suppose some mi are odd. Since m is even, all must be odd and 4 
does not divide m. Since all are odd, cos(mj7r/2) = 0 and the first product equals 1. 
The second product is obviously 0 and the bracket is again 0. If all are even, but 4 does 
not divide m, then exactly one of the mi/2 is odd. Therefore the first product is 0 and the 
second equals —1. Again the bracket is 0. Finally, suppose 4 divides m. Then 4 divides all 
mi, and we can reduce the problem to m' = m/4 and apply the same argument recursively. 

□ 


Remark. In dimension 5 we have 36 = 6^ = 4 x 3^. The first decomposition does give 0 in the 
bracket. However the second gives 1. 

In dimension 6 we have 36 = 6^ = 2 x 4^ + 4 x = 5^ + 2 x 2^ + 3 x = 5^ + 3^ + 2 x 1^. 
Hence in dimensions 6 and higher, any integer smaller than m may appear in the decomposition 
for m^. Therefore an argument based on divisibility will most likely fail. 


4.3. Nonnegative eigenfunctions. To prove Theorem 1.4 we will generalize the approach used 
on equilateral triangles in Section 3.2. We will show that px integrates to 0 over a union of small 
cubes with codimension one contained in one of the faces. Since an eigenfunction cannot equal 
0 on an open subset of the boundary (it already satisfies the Neumann condition there), it must 
change sign in the union of these cubes. Note also that it is irrelevant if these cubes are disjoint, 
but they must be subsets of the face. 

Suppose that \ = mf and consider the following set of points uniformly distributed on (—1,1). 

X = Ixk = ^— fc = 0, ...,m —ll 

( m J 


COs{l7lXk) 


k=0 


By [16, Section 1.341, Formula 3] 

'0 0 < / < m, 

m(—1)™+! I = rn. 

Now take a lattice of cubes with centers on and side length 2a. That is 
jC = = {y ■ yn = l,max \xi -yi\ < a} : X e 

Note that all cubes Cx are on one face of [—1,1]” if a < 1/m. 
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Consider one sequenee rrii and one permutation in the definition of ip\. The integral over the 
lattiee of the resulting funetion equals 



n—1 

cos(m<^(„)7r) n cos{mcr{i)'KZi)dzi... dzn-i = 

i=l 


n—1 m—1 

= cos(m^(„)7r) HE cos{m^^i)7rzi)dzi = 

i=l k=0 


n—1 m—1 

= cos(m<^(n)7r) HE 

i=l k=0 


-sin(m<^(i)7ra) cos{m^MnXk) = 


{ 0 , (mi,... ,m„) 7 ^ (0 ,... ,0,m), 

(-l)”"( 2 am)’""\ m„(ri) = m, 

i sin(m7ra)(—l)”^’''^(2am)"'“^, ma(i) = m for some i < n. 


( 8 ) 

(9) 


Note that in (8) we mean = 1 if a: = 0. The top ease in (9) is equivalent to /c > 0 in Seetion 3.2, 
while the other two eases eorrespond to the integral from (6). 

Henee 


Y1 / 

c^ec 



2ao,...,o,m(2am)” \n - 1)! am - 


n — 1 


■ sim 


71 


mTia) 


The last expression equals 0 if we ehoose 0<a<l/mso that 

Tram = (n — 1) sin(7r am). 

The existenee of sueh a is equivalent to the existenee of 0 < a: < vr sueh that 

X = {n — 1) sinx. 


( 10 ) 


This equation has a positive solution when n > 2. This proves that in dimensions n > 2 any eigen- 
funetion of a eube must ehange sign on the boundary. However this argument fails in dimension 
2, and Proposition 4.1 (or the earlier result [18] by Hoffman-Ostenhof) is the best we ean expeet. 
It is remarkable that (10) is exaetly the same as the equation for a the the equilateral ease (perhaps 
hinting at the faet that the equilateral triangle ean be embedded in a eube as an interseetion of that 
eube with a plane). 


4.4. General boxes. Consider an n-dimensional box with sides of length 2aj eentered at the ori¬ 
gin. The eigenvalues A for this box ean be indexed using a sequenee L of n natural numbers Zj sueh 
that 


A 


9 n 
7 T 

1 

2=1 




( 11 ) 


The eomplete set of eigenfunetions is given by 


n 

= YlF{li7iXi/ai), 

i=l 




11 



Figure 1. Regular heptagon deeomposed into subequilateral triangles, and reg¬ 
ular hexagon deeomposed into equilateral triangles. 

where F is either sine or eosine. However, any eigenfunetion that is nonnegative on the boundary 
ean be axially symmetrized by summing over all sign ehanges for all eoordinates. This proeedure 
still gives a nonnegative boundary and eliminates all oeeurrenees of sine. Therefore we ean assume 

that F{x) = cosx. 

Any eigenfunetion of a box, when restrieted to a faee, is also a sum of eigenfunetions on eaeh 
face (put Xi = tti for some i). This lower dimensional combination of eigenfunctions consists 
of eigenfunctions that are orthogonal to a constant eigenfunction (that is, they integrate to 0 over 
the face), and/or a constant term. If the constant term is not present, the linear combination must 
change sign on the face. Therefore, an eigenfunction that is nonnegative on the boundary must have 
a constant term when restricted to any face. Hence, its eigenvalue must admit indexing sequences 
Lj = {k = 6j{i)lj}. Taking L = Lj in Equation (11) thus yields 

4 af 4 a| 4 

This immediately proves that if any ratio of two squares of the side lengths is not the square of 
a rational number, then nonnegative eigenfunctions do not exist. 

For any / 7 ^ j we have 


Hence ai/aj is also rational for any i ^ j. Therefore ai = r^a for some rational r* and real a, and 
this box can be used to tile a cube. Then any eigenfunction positive on the boundary of this box 
gives an eigenfunction on a cube with the same property (thanks to Neumann boundary matching 
in the tiling). But we proved these do not exist. Therefore Theorem 1.4 also holds for arbitrary 
boxes. 


5. Proof of Theorem 1.3 for n > 6 

For a regular hexagon we can simply take the symmetric mode of the equilateral triangle 
(defined in (3)) and cover the hexagon with its reflections to get an eigenfunction which is positive 
on the boundary. 
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Figure 2. Regular pentagon decomposed into acute superequilateral triangles, 
blue triangle T and red rhombus R. 

Now consider a regular polygon with n sides, where n > 6. Such a polygon can be decomposed 
into n subequilateral triangles {ABD on Figure 1). The second Neumann eigenvalue of a subequi- 
lateral triangle is simple (Lemma 2.1) and the second Neumann eigenfunction is symmetric [24, 
Theorem 3.1]. Hence it is also the second eigenfunction on the right triangle formed by cutting the 
isosceles triangle in half {ABO and ADO on Figure 1). 

The second Neumann eigenfunction must have exactly 2 nodal domains, by Courant’s nodal 
domain theorem (see e.g. [13, Sec. V.5, VI.6]). By symmetry, the nodal line must either connect 
the two long sides {AB and AD) of the subequilateral triangle, or start and end on the short side 
(BD). From Lemma 2.3, the second case is not possible, regardless of the shape of the triangle. 
Hence this eigenfunction is positive on the short side, and it can be reflected n times inside of the 
regular polygon to cover the whole regular polygon. We obtain an eigenfunction on the regular 
polygon that is positive on the boundary. Therefore Theorem 1.3 holds for n > 6. 

As a corollary from the above proof we also get a partial result for Conjecture 2.4 

Corollary 5.1. The nodal line for the second Neumann eigenfunction on right triangles with small¬ 
est angle a < ti/Q connects the interiors of the two longest sides. 

6. Proof of Theorem 1.3 for regular pentagons. 

A regular pentagon decomposes into acute superequilateral triangles instead of subequilateral 
triangles (as was the case of n > 6 sides). 

The second Neumann eigenvalue /i 2 of a superequilateral triangle {ABD on Figure 2) is simple 
but the second eigenfunction is antisymmetric [24, Theorem 3.2] (as opposed to symmetric for 
subequilateral triangles). By Lemma 2.2 all eigenfunctions for /is are therefore symmetric. But all 
these eigenfunctions belong to the second (simple) eigenvalue of the right triangle OAB obtained 
by cutting the isosceles triangle ABD in half (shaded on Figure 2). Therefore /is of a superequilat¬ 
eral triangle ABD is simple, with the eigenfunction symmetric with respect to OA. Unfortunately, 
Lemma 2.3 applies only to eigenfunctions for /i 2 . Moreover, we need to exclude a possibility of 
having 3 nodal domains (allowed for /is). 

These two problems make the pentagonal case much harder than regular polygons with n > 6 
sides. Moreover, there is essentially no hope of finding explicit trigonometric formulae for its 
eigenfunctions. A completely different approach is required. 
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Consider the rhombus R {ABCD on Figure 2) built using right triangle T {ABD on the same 
figure) with the smallest angle at least 7 r /6 (equal to vr/S for our regular pentagon). Then [34, 
Corollary 1.3] gives 

fj^iiR) < Ai(i?). 

Note that the elassieal Levine-Weinberger inequality [25] only gives /is < Ai. Furthermore, the 
eigenfunetion U 2 that belongs to jJi 2 {T) extends by symmetry to a doubly symmetrie eigenfunetion 
u on R. Then u must belong to the lowest eigenvalue of R whieh possesses a doubly symmetric 
mode. Otherwise, a doubly symmetric eigenfunction of the lower eigenvalue of R would be an 
eigenfunction for T. Therefore 

fi2{T) = fJ,4{R)- 

For the triangle T = ABD we have 

Lemma 6.1. The partial derivatives and Uy of the second Neumann eigenfunction u of T are 
never zero and have opposite signs. 

Remark. Note that this result can be deduced from the last paragraph on page 244 of Atar- 
Burdzy [6]. Nevertheless we present a simpler proof. 

Proof. We will follow the proofs of [28, Lemmas 3.2,3.4] and [33, Theorem 2]. First note that [33, 
Lemma 2] applies to T, hence its second Neumann eigenfunction u is strictly monotonic on AB. 
We can assume that Ux> Q and < 0 on AB. 

Now we consider the doubly symmetric extension of u to the rhombus R. On CB we also have 
Uy < due to double symmetry of u, while on CD and DA we have Uy > 0. Similarly, Q 
on DA, and on CB and CD. Furthermore, is antisymmetric with respect to r/-axis 

and symmetric with respect to x-axis (again by double symmetry of u), while Uy has reversed 
symmetries. 

Suppose Uy is zero somewhere in R, then by [28, Proposition 2.1(i)] it must change sign inside 
R. By antisymmetry, it must be positive somewhere in ABC. But My < 0 on AB and CB and 
My = 0 on AC. Hence a nodal domain of My is a subset of ABC (part of the nodal line might be a 
subset of AC). We already noticed that u belongs to ^^{R), hence 

PiiR) = Ai(A^) > \i{R) > PiiR), 

a contradiction. Hence My < 0 on ABC (hence also on T). Similarly we can prove that m^ > 0 on 
T. □ 

We need a domain monotonicity result for the eigenvalues of the domains with mixed boundary 
conditions. This is a special case of a more general partial domain monotonicity principle proved 
by Harrell. We present this special case due to its rather simple proof. 

Lemma 6.2 (Special case of Harrell [17, Corollary 11.2]). Suppose Di C D 2 are open and the Neu¬ 
mann boundary d^Di of Di is contained in the Neumann boundary di\]D 2 of D 2 (see Figure 3). 
Then the lowest eigenvalue on D 2 for the mixed Dirichlet-Neumann problem is smaller than the 
lowest mixed eigenvalue on Di, unless Di = D 2 and d^Di = djqD 2 . 

Proof. Suppose p is the eigenfunction for Di. Extend it with 0 to the whole set D 2 . Note that 
the extension satisfies the Dirichlet boundary condition on dD 2 \ OnDi, hence on the Dirichlet 
boundary dDD 2 . Note also that d^Di does not intersect D 2 \Di, hence the extension is continuous. 
Therefore it is a valid trial function for the Rayleigh-Ritz quotient on D 2 . But it also equals 0 on an 


14 


OnDi 



dDD2 


Figure 3. Domain monotonicity: Di c D 2 and Neumann boundary condition 
on Di is specified only on a portion of the Neumann boundary 0^02- Solid lines 
indicate a Neumann boundary while dashed lines indicate Dirichlet boundarh. 

open set (if Di strictly included in D 2 ), or equals 0 on a piece of boundary 8^02 \ d^^Di (satisfies 
both Dirichlet and Neumann condition). In either case it must be 0 everywhere. □ 

Corollary 6.3. Let Di C D be a nodal domain for the eigenfunction for fi 2 {D). Suppose we find 
D 2 <Z D such that d^Di C Oat £*2 C dD, and the mixed eigenvalue \ fiD 2 ) > P 2 {D). Then the 
nodal line for p 2 {D) must intersect the Dirichlet boundary c)d-D 2 - 

Proof The mixed eigenvalue of Di equals p 2 {D) and is smaller than \ fiD 2 ). Note also that for 
a nodal domain we always have d^Di C dD and doDi is the nodal line. If Di C D 2 , then the 
above lemma gives Ai(Zi)i) > A 1(^2), leading to a contradiction. Hence Di (f D2, and the nodal 
line OdDi must have a nonempty intersection with D \ D 2 . Hence it must intersect c)d-D 2 - D 

Lemma 6.4. The lowest eigenvalue of a mixed Dirichlet-Neumann eigenvalue problem is simple 
and the eigenfunction can be taken positive inside of the domain. 

Proof. Let u be the eigenfunction for the smallest mixed eigenvalue. Then |m| gives the same 
value of the Rayleigh-Ritz quotient. This means that |m| is also a minimizer of the Rayleigh-Ritz 
quotient. Any minimizer is an eigenfunction (see e.g. [23, page 55]). Then elliptic regularity shows 
that |m| solves the eigenvalue problem classically (not in the weak sense). But A|m| = —A|m| < 0, 
so |m| is superharmonic. Take any point p inside the domain at which |m| = 0. Nonnegativity at p 
violates the minimum principle. Therefore m > 0 inside of the domain. 

Suppose u > 0 and v are orthogonal eigenfunctions belonging to the lowest eigenvalue. Orthog¬ 
onality forces V to change sign in the domain, and the argument above shows that |r;| is yet another 
eigenfunction, violating the minimum principle. Therefore the eigenvalue is simple. □ 

Now we restrict our attention to the right triangle T = OAB with angle vr/S near A (see 
Figure 2). We can assume that \OA\ = 1 and the second Neumann eigenfunction u is negative 
at B. We know that is positive and Uy is negative by Lemma 6.1. This implies: 

Lemma 6.5. The nodal line for the eigenfunction for the lowest positive eigenvalue on the right 
triangle OAB is the graph of a strictly increasing function, hence it must touch the longest side 
AB. 

Suppose that the nodal line also touches the shortest side OB (or hits the origin O). We construct 
an explicit subdomain D 2 of OAB, such that a nodal domain of the second eigenfunction of T is 
contained in D 2 , yet Ai(iA 2 ) > P 2 {T), contradicting Corollary 6.3. 

Therefore the second eigenfunction of T has a fixed sign on OB. Reflecting this eigenfunc¬ 
tion 10 times we can cover a regular pentagon and find an eigenfunction which is negative on its 
boundary. This proves the remaining case n = 5 in Theorem 1.3. 
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Number of nodes N 

Rayleigh Quotient 

38 

12.2482 

128 

12.2476 

463 

12.2475 


Table 1. Upper bounds for fiT' using conforming quadratic finite elements (vali¬ 
dated numerics. 


Lemma 6 . 6 . Let fl be a Lipschitz domain and fit = {(a^, ty) ■ {x, y) E fl} (a stretched domain), 
fort > 1. Then for any mixed boundary conditions \i{flt) < Ai(U). The same is true for the 


smallest nonzero Neumann eigenvalues. 

Proof See the last paragraph on page 132 of [24]. □ 

To avoid irrational triangle vertices, we define right triangles T' C T C T", as well as a 
reference triangle Tref. These triangles have two common vertices O = (0,0) and A = (1,0), 
and their third vertices are as follows: 

T' :B' = (0, 85/117) = (0, 0.72649573 ...), (12) 

T : S = (0, tan(7r/5) = (0, 0.72654252 ...), (13) 

T" : B" = (0, 93/128) = (0, 0.7265625), (14) 

T/ef :5ref = (0,1). (15) 

The triangle T can be used to create the regular pentagon. Let us define linear maps 

T • ^ref T, (16) 

f : Tref —^ T" (17) 

which fix A and send Tref to B and B", respsectively. 

Lemma 6.7. The smallest nonzero Neumann eigenvalues ofT and T' satisfy 

Pt ^ pT' ^ 12.25. 


Proof The triangle T is obtained by vertically stretching T', hence Lemma 6.2 gives the required 
monotonicity. We can get a very accurate upper bound for pt' using the finite element method with 
quadratic nodal conforming elements. An eigenfunction approximation obtained from conforming 
elements is a valid trial function (continuous, piecewise quadratic) for the Rayleigh-Ritz quotient, 
hence we get a strict upper bound for px' by plugging it directly into the Rayleigh-Ritz quotient. 
We tesselate T' by congruent triangles with N distinct nodes, and construct test functions in (T') 
which are continuous and piecewise quadratic. Standard finite element approximation arguments 
tell us that these Rayleigh-Ritz quotients will form a decreasing sequence in N, which converges to 
the true p(T'). In Table 1 we show the Rayleigh-Ritz quotient computed using interval arithmetic 
(we only report 4 digits, though the interval around these approximations is of width le — 9). We 
can comfortably bound p{T') above: 


Pt' < 12.2483 < 12.25. 


□ 
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In what follows we will use a threshold value V = 12.25. We choose this particular value to 
make sure <V- 

Lemma 6.8. Let D C T^^/. Clearly'ip{D) C T",ip{D) C T. Impose any mixedDirichlet-Neumann 
boundary conditions on dD. Then ijj and p induce mixed boundary conditions on 'ip{dD) and 
(p{dD). The smallest eigenvalues for the mixed eigenvalue problems on <p{D) and 4>{D) then 
satisfy 

\{p{D)) > X{f{D)). 

Proof The domain f{D) is a vertically stretched version p{D). Therefore Lemma 6.6 implies the 
result. □ 

Definition (cf. Figure 4). Let a = 1 and b = tan(7r/5) and N = 64. Define a grid on T: 

G = {{ai/NM/N) 

and for a grid point p let 

Rl{p) = {ix,y) \x>xo,y< yo} fl T, 

Ru{p) = {ix,y) \ x <xo,y> yo} fl T. 

Furthermore define two sequences of points: 

pz.= {(28,21), (31,27), (24,14), (28, 23), (26,19), 

(21,9), (30,27), (23,15), (31,30)}, 

{(25,24), (20,17), (28, 30), (18,13), (22,18), (28,28), 

(29,29), (19,12), (27,25), (25,20)}. 

ik) (k) 

For these particular points p\ and p\j we define sequences of subdomains of T. 

Definition. Let = 0, and 

L{k+i) =L'=Ui?L(pi?'+^^), 

We call these sets the upper and lower exclusion regions respectively, for reasons which will be¬ 
come clear shortly. Also define the following subdomains ofT with associated mixed boundary 
conditions: 

Du{k,p) = T \ (L*-^^ U Rl{p)), Dirichlet on dDjj{k,p) \ {AB U OB),Neumann elsewhere 

( 20 ) 

Dcik^p) = T\ U Ru{p)), Dirichlet on dDL{k,p) \ dT, Neumann elsewhere (21) 

(fc + l) (fc + l} 

Remark. We note that for appropriate choices of points p = p} \Pu Equation 18 and 
Equation 19 will make 

=T\(L('=+')), Zlz.(fc,pS?+'^) =T\(f/('^+i)). 


(18) 

(19) 


We use will this fact later on. 
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Lemma 6.9. The lowest mixed eigenvalues on the following subdomains ofT and T” = fxp 
satisfy 

Xi(^Du > V, for0<k<9, (22) 

Ai(Di(fc + l,pr'^)) >A", for0<k<8, (23) 

Ai (Di(10, (0,0))) > Ai {fp-^DLilO, (0,0))) > )/. (24) 

Remark. Note that in the last case we use p = (0, 0), which results in no enlargement to 
Effectively, ^^(10, (0, 0)) = T \ 

The proof is postponed to Seetion 8. The proof involves lower bounds for the eigenvalues of 
large sparse rational matriees. 

Lemma 6.10. Assuming the nodal line for pT starts on the side OB, it does not intersect any of 
the excluded sets and (e.g. the shaded regions on Figure 4). 

Proof. The elaim is elearly true for and L'^^X Assume the elaim holds for and L^^X 
We eonstruet Djj{k,p^^~^^^) and note that Xi{Du{k,p^^^^^)) is larger than the threshold V (by 
Lemma 6.9). Henee L)u(k ,is a subdomain of T whose lowest eigenvalue exeeeds the 
smallest nonzero Neumann eigenvalue px for T. Therefore Corollary 6.3 implies that the nodal 
line for the eigenvalue px must interseet Du{k,p^~^^'‘). This cannot happen on L^’^X and hence 
the nodal line intersects Rx{p^'^^^) \ L^^X Monotonicity of the nodal line (Lemma 6.5) ensures 
that point p^x^^^ separated from L'X) by the nodal line. In particular it must belong to a different 
nodal domain than as does any point above and to the left of p^x^^\ Hence the nodal line 
cannot intersect the constructed U^^^^X 
A similar argument applies to L^^^^X 

□ 


Theorem 6.11. The nodal line for the eigenfunction for px does not intersect OB. 

Proof. Suppose the nodal line has its endpoints on AB and OB (Lemma 6.5 ensures it ends on 
AB). Divide T into two nodal domains of the eigenfunction of px- By construction, the smallest 
mixed eigenvalues of these nodal domains equal px, which is not larger than the threshold V. 
Furthermore,one of these nodal domains is contained in T \ U^^^X 

Domain monotonicity (Lemma 6.2) and the last case of Lemma 6.9 lead to a contradiction. □ 

This theorem ensures that the eigenfunction for px of the right triangle T does not change sign 
on the shortest side. Its symmetric extension to the regular pentagon is an eigenfunction of that 
pentagon, which does not change sign on the boundary, proving the remaining case of Theorem 1.3. 

7. The exclusion search algorithm 

Our proof for pentagons contains a rather strange collection of domains and associated grid 
points Px and pu defined in (18) and (19). In Lemma 6.9 we check that the smallest mixed eigen¬ 
values of these domains satisfy appropriate inequalities. The lemma does not depend on, or explain, 
how these exact domains were chosen. In reality, these domains were found by a computer search. 
In this section, we present an algorithm which generates these points px and pu and the associated 
domains based on numerical computations. In the following section we describe how, for a given 
domain, we validate the assertion of Lemma 6.9. 
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The proof in the preeeding seetion is a human-readable, computer-assisted proof, and it is worth 
making a comment about the structure of such proofs. Of course, all mathematicians appreciate 
a short, beautiful proof of a theorem, but it is not clear what standard of “beauty” should apply 
to a computer’s work. In this section, we argue that a computer proof should strive to have small 
certificates, be easily checkable, and be adaptable to other situations, and that our proof has these 
properties. We direct the interested reader to [29, Ch. 2] for a reference on human-readable, 
computer-assisted proofs. 

7.1. Small Certificate. The computer-generated part of a human-readable proof is called the 
proof certificate. Because computers can generate a lot of output very quickly, the proof cer¬ 
tificates that they generate can be very long and time-consuming to read, unless care is taken to 
program the computer to search for short proofs. 

For example, Appel and Haken’s famous proof of the four-color theorem in graph theory [1, 2, 3] 
is partly computer-generated; its certificate consists of a certain argument by cases. This portion of 
the argument is very long: it occuped a lengthy microfiche supplement to the articles [1,2]. A later 
version of the proof appeared in the book [3], which is over 700 pages long. A mathematician who 
wanted to read their proof would, amongst other tasks, have to inspect all of this work carefully; it 
is likely that this would take years of full-time work to do. More recent research on the four-color 
theorem, such as [30], has in large part focused on finding computer-generated proofs with shorter, 
simpler proof certificates. These proofs are much less time-consuming for a mathematician to read. 

We found a rather short proof certificate for our problem: it is the two lists of points (18), (19), 
together with the upper bound px' < 12.2483 in the proof of Lemma 6.7. 

7.2. Checkability. When possible, a computer-assisted proof should be checkable without requir¬ 
ing specialized software. To check our proof, a human needs only use standard, validated numerics 
sofware (such as I NT LAB [3 1 ], for interval arithmetic) to give bounds on the eigenvalues of certain 
matrices which we describe. An interested reader can use this, or other standard interval arithmetic 
software, to check the proof without needing to develop their own specialized code. 

7.3. Adaptability. A mathematician should be able to adapt a computer-assisted proof to other 
problems, without requiring significant changes in the human-written part. Our proof has this 
advantage: it is possible to change the grid size N = 64, or triangle height b = tan(7r/5), and 
re-run our algorithm. This would produce a new list of points pi and pu, and a new upper bound 
for T'. Having made these changes, the reader would select a new threshold V exceeding this 
upper bound, and then verify the assertions of Lemma 6.9 in order to establish the analogue of 
Theorem 6.11. 

7.4. Tbe algoritbm. We now present a strategy. Algorithm 1, for constructing certain subdomains 

Dl, Du, fc = 0,.. of T. These are the subdomains for which we shall study the low¬ 

est mixed eigenvalues. The algorithm fails if both of the subdomains and cannot be 
constructed for some k. The steps of the algorithm are shown on Figure 4. 
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(a) Steps (1) and (2): shaded with red, 

magenta and Ru green. Dirichlet boundary condi¬ 
tions on domain Du{k,p) as red lines, Neumann 
as blue lines. Test domain DL{k,p) is the com¬ 
plement of magenta and green. 


(b) Steps (3) and (4): shaded with red, 

magenta and Rl green. Dirichlet boundary 
conditions on domain D^^k -f l,p) as red lines, 
Neumann as blue lines. Test domain Du{k,p) is 
the complement of magenta and green. 


Figure 4. Algorithm for constructing excluded domains. 

Algorithm 1 Construction of of subdomains U and L. 

1: Begin with and and set threshold value V as in the previous section and let A: = 0. 

2: for k=0,l,2.. do 
3; Let S’ = 0. 

4; for all points p on the grid G (See Remark 1 below) do 
5; Construct the domain Df/(A:,p). 

6: if Xi{Du{k,p)) > V then, add p to the set S. 

7; if S is empty then, := = (0,0) 

8 : else 

9; for all points p G S' do construct the domain DL{k,p). 

10: Let be the point p such that the smallest mixed eigenvalue of Di{k,p) is maximal (over 

S). See Figure 4a. Point defines . 

11: Let 5 = 0. 

12: for all points p G G (See Remark 1 below) do 

13: Construct the domain DL(k + l,p) (we are using the newly created 

14: if Ai(Z9l(A; + l,p)) > V then add p to the set S. 

15: if S is empty and it was empty at step 7 then return ’Algorithm failed’ and exit 

16: if 5 is empty then, := = (0,0) 

17: else 

18: for all points p G 5 do construct the domain Du{k,p). Here we are using 

19: Let be the point p such that the smallest eigenvalue of Djj{k,p) is maximal. See 

Figure 4b. Point defines 

20: if The smallest mixed eigenvalue for or Djj{k,p^^~^^^) exceeds V then 

21: return ’Algorithm is successful’ and terminate. 

22: else A: •(— A; + 1 and go back to step 2. 
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Remarks. 

(1) Note that it is not necessary to construct Du{k, p) for every grid point p. Indeed, ifp G S 
then any point from Ru{p) is also in S, by domain monotonicity Lemma 6.2. Similarly if 
p ^ S, then neither are any points from Riip)- As a consequence, we only perform a linear 
search starting at the origin and moving up until we find the first grid point p E S, then 
move one grid-point right and repeat. 

(2) Steps 2 and 4 choose the “best” pointp to use to expand the excluded regions, with the aim 
of fine-tuning the algorithm to terminate in as few steps as possible. It is entirely possible 
that the length of the lists in (18) and (19) could be reduced by using different strategy. On 
the other hand, the above formulation is simple, and more elaborate choices we tried did 
not lead to significant improvements. 

(3) We choose to grow the exclusion sets U and L by adding only one point p. We could simply 
add all points p E S to new exclusion sets. However, for every p we use, we need to supply 
either an exact eigenvalue or a validated lower bound. In other words the lists of points 
Pl and pu could be replaced with lists of lists of points, and in every step the whole list 
of points would be used to grow and This is certainly better from a numerical 
point of view, but this improvement would come at the expense of an extremely long proof 
certificate. 


8. Validated lower bounds 

There are two teehnieal diffieulties, whieh we resolve here, in obtaining rigorous numerieal 
bounds for eigenvalues. By eonsidering T' C T C T" (see (12-14)) with rational eoordinates in 
T' and T" we avoid floating point errors due to inexaet domain representation. We further map 
T" —)■ T"', a similar triangle with perpendieular sides 93 and 128. If A is an eigenvalue of T"', then 
128^A is an eigenvalue of T". 

Before we deseribe the validated numerieal strategy, we list the sub-domains of T"' for whieh we 
will eompute the lowest mixed eigenvalues. Next, in Seetion 8.2, we fix notation and definitions 
for the numerieal approximation strategy used. 

8.1. Polygonal domains for which we find validated lower hounds. We now enumerate the 
upper and lower domains, by listing their vertiees. Note that for ease of reading these are listed on 
the (64, 64) right isoseeles triangle T^e/. To obtain the eoordinates of these points in T"', seale the 
x-eoordinate by 2 and the y-eoordinate by 2 q; = 93/64. 

The upper domains Djj{k, (shown on Figure 5) are created using grid points pjj. We only 

list every other vertex on the step-like (red) portion of the boundary of the polygon and we skip the 
top vertex (0, 64) (connecting two blue boundary pieces): 

A: = 0 : (0, 0), (25, 24) 

A: = 1 : (0, 0), (20, 17), (28, 21) 
k = 2 : (0, 0), (28, 30) 

A: = 3 : (0, 0), (18, 13), (24, 14), (28, 21), (31, 27) 

A: = 4 : (0, 0), (22, 18), (28, 23), (31, 27), 

A: = 5 : (0, 0), (24, 14), (26, 19), (28, 28), 

A: = 6 : (0, 0), (21, 9), (24, 14), (26, 19), (28, 23), (29, 29) 

A: = 7 : (0, 0), (19, 12), (24, 14), (26, 19), (28, 23), (30, 27) 

A: = 8 : (0, 0), (21, 9), (23, 15), (26, 19), (27, 25), (30, 27) 

A: = 9 : (0, 0), (21, 9), (23, 15), (25, 20), (28, 23), (30, 27), (31, 30) 
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Figure 5. Upper domains from the algorithm. 



Figure 6. Lower domains from the algorithm. 


Lower domains DL{k + (shown on Figure 6) are ereated using points pi- We again 

skip some vertiees whieh ean be deduees from the rest. 

k = 0 : (28, 21) 

k = l : (20, 17), (25, 24), (31, 27) 

k = 2 : (24, 14), (25, 24), (28, 30) 

k = 3 : (18, 13), (20, 17), (28, 23) 

k = A : (18, 13), (20, 17), (22, 18), (26, 19), (28, 30) 

k = 5 : (21, 9), (22, 18), (25, 24), (28, 28) 

k = 6 : (18, 13), (20, 17), (22, 18), (25, 24), (30, 27) 

k = 7 : (19, 12), (23, 15), (25, 24), (28, 28), (29, 29) 

k = 8 : (19, 12), (20, 17), (22, 18), (25, 24), (27, 25), (28, 28), (29, 29), (31, 30) 

k = 9 : (19, 12), (20, 17), (22, 18), (25, 20), (27, 25), (28, 28), (29, 29) 
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8.2. Nonconforming finite elements method. Suppose is any of the polygonal domains listed 
above in Section 8.1. Let Fd be the Dirichlet part of its boundary. Since there is no closed- 
form expression for the eigenvalues, a lower bound for the smallest mixed eigenvalue Ai must be 
determined with the aid of approximation theory. 

We can reformulate the mixed eigenvalue problem for {uj, \j) on the domain in variational form: 
find the jth eigenpair {uj, \j) so that for all test functions v in the Sobolev space 


Vuj ■ Vu = Xj 


UjV. 


Choose a suitable finite-dimensional function space 14 (where > 0 is called the mesh parameter), 
and consider an approximation Uj^h G 14 of Uj, which satisfies the discrete eigenvalue problem-. 
find the jth discrete eigenpair (wyh, Xj^h) so that for all test functions Vh G 14, 


/ ^j,h I '^j,h'^h- 

Jn Jn 

By choosing Vh to be basis functions of 14, we obtain a discrete generalized eigenvalue system of 
the form 

^h^j,h Aj/jM/jUj/j (25) 

The matrix Kh is called the stiffness matrix and the matrix M/i is the mass matrix. 

Recall that if the integrals in (25) are computed exactly, and 14 C iFp^(f2), then the method is 
called conforming, otherwise it is non-conforming. For a conforming method it is easy to see from 
the definition of the Rayleigh-Ritz quotient that Ai < Ai/j. However, we want a lower bound for 
Ai. For this, we shall use a non-conforming method (following [4, 11]). In what follows, we will be 
using a the well-known Crouzeix-Raviart linear non-conforming finite element discretization. Let 
T = QCD denote a right-angled triangle which is similar to OAB", with a right angle at Q. Let 
\QC\ = h, \QC\ = a/iwherea = In what follows h = 2. We tesselate a polygonal subdomain 
H C OAB" by a mesh T = T{h) consisting of congruent copies T of T, ie, OAB" = 

Recall that a linear function on T can be uniquely described by specifying its values at three points 
(these are its ’degrees of freedom’). For the Crouzeix-Raviart method we use, these degrees of 
freedom are on the (rational) midpoints of the edges of T. Let 8 denote the collection of all edges 
in the tesselation, and let F d denote the Dirichlet boundary of D. The finite-dimensional function 
space we use is 

Vn = CR\){T) -.= {u e L^(D) 44 e Pi(T), u is continuous at interior nodes 
and 0 at the nodes on £^(Fp))} 


It is obvious, now, that with this choice of 14 the matrices in (25) are symmetric and that all 
entries are rational. Also, the mass matrix is a diagonal matrix with diagonal blocks ^Isxs- 
In our tesselation, it is clear that each node Xk is adjacent to at most 4 other nodes in T (recall that 
in the Crouzeix-Raviart discretization, each interior node lies at the midpoints of the sides of two 
triangles in the mesh, and its adjacent nodes are the midpoints of the other sides of those triangles). 
Therefore, K/j will have atmost 5 non-zero entries per row and column. By examining the local 
stiffness matrix for T G fl, 

2/a —2/a 0 \ 

—2/a 2(Q;-|-l/a) —2a 1 . 

0 —2a 2a J 


[Kh] = 


Vfi ■ V(f)j 


Jt 
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The entries of the full stiffness matrix are built from these loeal eontributions. The diagonal 
entries of K/j are, therefore, equal to 2a, 4a, 2/a, 4/a, 2(a + l/a), or 4(a + l/a). The off-diagonal 
terms are equal to ±2a or ±2/a. The simple strueture of allows us to write 




h'^a 


^j,h^h^j,h 


(26) 


where Bh is diagonal with 1 or 2 as the only nonzero entries. This matrix Bh is easy to faetor 
and to invert. We also define a sealed version of the stiffness matrix whieh has integer entries 
Ch ;= (93)(64)Kh. We ean then rewrite (26) as 




(93)(64) 


h'^a 


93\ 


(27) 


We denote by Bh := (B/i) This is a diagonal matrix with 1 or 1/\/2 on the diagonals. With 
this notation we ean rewrite the generalized eigenvalue problem in (26) to 

hj^oi 

^h^j,h — (b/iK/jD/j)Uj/j ^ Xj fiMj fi = Ij iiMj ii (28) 

The scaling in h above is typically important from the point of view of performing stable compu¬ 
tations; the conditioning of the original discrete system deteriorates with h. However, note that we 
use the fixed value h = 2. 

Each of the stiffness matrices Kh and mass matrices for T'" and for the various subdomains 
Du(-), Dl{-) has roughly 4000 rows and columns, so we do not reproduce these matrices here. 
Rather, we describe how the matrices were constructed. We first used FEniCS [26] to assemble the 
stiffness and mass matrices on T”' (the full triangle, not the subdomains). We rescaled the resulting 
matrices to get integer entries, and transferred the results to a scipy. sparse [20] matrix for 
further processing, including handling of boundary conditions. For the subdomains, to enforce 
the shape of a subdomain and the Dirichlet condition, we set the relevant degrees of freedom to 
0, forcing the solutions to equal 0 on the complements of the subdomains from Lemma 6.9. In 
practice, we accomplish this by removing the rows and columns for these degrees of freedom from 
both mass and stiffness matrices. 


8.3. Generalization of the validated nonconforming lower bounds. We now describe how to 
transform the above numerical scheme into rigorous lower bounds for mixed Dirichlet-Neumann 
eigenvalues. First, for the purpose of validated numerics, the entries of the matrices K^, should 
be rational numbers. This can be done by ensuring the vertices of the domains above are rational, 
and by defining 14 appropriately, which was done in the preceding section. Next, the eigenvalues 
of large discrete systems such as (25) cannot be found in closed form. We must therefore ensure 
any approximation errors incurred are controlled in a manner that the desired lower bound is robust. 
We present three approaches in Section 8.4 and Section 8.5 below. 

We begin by surveying the literature on such rigorous lower bounds. In [4], an asymptotic result 
shows that provided h is small enough, Xi^h computed using a (specific) non-conforming method 
provides a lower bound for the first eigenvalue Ai for the Dirichlet problem. It is only recently 
that a computable lower bound based on a non-conforming finite element method for the first 
Dirichlet eigenvalue of the Laplacian on a polygonal domain has become available, in [11]. This 
is achieved by approximating the first Dirichlet eigenfunctions using piecewise linear Crouzeix- 
Raviart elements, that is, approximating the first eigenfunction by functions from 

CRq{T) := G LF‘{Vt)\v\T G Pi(T), n continuous at the midpoints of interior edges 

and 0 at the midpoints of boundary edges} 
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Suppose CRq{T) is used to construet the generalized eigenvalue problem of form (25) for the pure 
Dirichlet problem. Barring very speeial eireumstances, the eigenvalues of this system cannot be 
found in closed form. Therefore, suppose we use an iterative approximation technique to obtain 
an approximation \cR,i to \i^h (which is, itself, an approximation to Ai). The authors of [11] 
establish the computable inequality Theorem 3.1 of their paper: if {\cr,1iUcr,i) G M x CRq{T) 
approximate the first eigenpair (Ai, Mi) (with eigenfunction normalized to 1), if r is the algebraic 
residual, and if Xcr,i is closer to the first true discrete eigenvalue Xcr,i than to any other discrete 
eigenvalue (see [11, Lemma 3.8]), then 


-^-< Ai. (29) 

1 + R?{XcR,i - |r|B-i)if2 

Here, H is the maximal diameter in the regular triangulation T and k is a universal constant 
given by = I + < 0.1931 for the first positive root of the Bessel function of the first 

kind. 

In our problem, we need to establish a similar lower bound for the first eigenvalue of the Lapla- 
cian with mixed data. We need to therefore extend Theorem 3.1 from [11] to this case. Specifically, 
we seek approximations ucr,i G CR]y{T) for the first eigenfunction in a mixed eigenvalue prob¬ 
lem. 

The non-conforming interpolant used to establish the interpolation estimate of Theorem 2.1 in 
[11] remains unaltered except for the obvious change in domain and range: X^c ■ CR^{T) 

is defined on all edges E which are not on T^ as: 


lNcv{mid{E)) 


1 

¥ 



V ds. 


Theorem 2.1 also holds, with n unaltered. Recall that this n is the optimal Poincare constant on a 
triangle [24]. Provided T^) is non-empty, the H^{D) semi-norm remains a full norm (a fact used 
in the argument of [11]). Therefore, the same inequality (29) holds. 

Note, the lower bounds by Carstensen and Gedicke are conditional on algebraic eigenvalue 
approximation being close the smallest true eigenvalue (see [11, Lemma 3.8]), and this condition 
is not easy to validate. Nevertheless, their Theorem 3.2, relating exact discrete eigenvalue Xcr,i 
and Ai, generalizes unconditionally to mixed eigenvalues: 


Lemma 8.1. Let (Acr,i) be the true eigenvalue corresponding to a shape-regular Crouzeix-Raviart 
discretization of the mixed eigenvalue problem for the Laplacian on a polygonal domain, with 
maximal mesh size H. Let Ai be the first eigenvalue of the mixed eigenvalue problem on this 
domain. Then 


XcR,l 

1 + R^XcRyH^ 


< Ai. 


(30) 


Moreover, if Xcr,i is an approximation to the discrete eigenvalue Xcr,i which is closer to Xcr,i 
than to XcR, 2 , then again 


_ XcR,l — klg-i_ 

1 -f k‘^{Xcr,i - 


< Ai. 


(31) 
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Furthermore, the left hand side of (30) is inereasing in A allowing us to find a lower bound for 
Ai using a lower bound for the diserete eigenvalue Xcr,i- In partieular 


A < XcR,l 


A 

TTr^aip 


< Ai. 


(32) 


Note that we have = {h/GA^ + {ah/GA^ and < 0.1931 (reealling that we take h = 2 and 
a = 93/128). To obtain a proof of Lemma 6.9, we need V = 12.25 < Ai. Therefore it is enough 
to show that Acj?,i > 12.25 + 3/256 =: A. This ehoiee of A is made sueh that if we multiply both 
sides of the diserete eigenvalue problem by a power of 2, then all matriees involved in (26) have 
integer entries. No floating point errors are ineurred at the stage of matrix assembly. 

As is standard for diserete eigenvalue eomputations, the quality of approximation to the first 
diserete eigenvalue Xcr,i depends (inversely) on the speetral gap between the first and seeond 
diserete eigenvalues, Xcr,2 — Xcr,i- 


8.4. Lower bound for the smallest eigenvalue of a sparse matrix. Here we present a method 
of finding a verifiable lower bound for the lowest eigenvalue of a sparse rational matrix, without 
eomputing any eigenvalues. We use this method on the stiffness and mass matrix systems deseribed 
above, and the domains from Lemma 6.9. 

We have several generalized eigenvalue problems Kx = Xcr,iMx, one for eaeh domain from 
Lemma 6.9. We want to show that the smallest eigenvalue of eaeh exeeeds A. Instead of finding 
an approximate lowest eigenvalue (and eheeking that it is aetually the smallest), we may prove that 
K — AM is positive definite. K — AM is positive definite if and only if \/M~^KVM~^ — AI is 
positive definite. If the latter is positive definite than the smallest eigenvalue of '/M^K'/M^ is 
larger than A, and it is the same as the smallest eigenvalue for the generalized system Kx — 'jMx. 

Proposition 8.2. The matrices K — AM are all positive definite. 

Proof. Sinee A has been ehosen to be a rational number, this is a routine eheek with standard soft¬ 
ware. For the diligent reader who wishes to eheek this on their own, there are several teehniques, 
all of whieh we have tried sueeessfully: one ean use the LU and LDL' deeomposition in interval 
arithmetie, or one ean eompute an exaet LU deeomposition with rational arithmetie. Here are the 
details of how to eheek the proposition with standard software. 

The matrix K — AM has rational entries and is sparse. We used CHOLMOD [12] to find a fill- 
redueing permutation (a matrix P whieh makes the LDL' deeomposition of P{K — 'yM)P' as 
sparse as possible). Then we find the LDL' or LU deeomposition of the permuted matrix. A 
matrix is positive definite if all diagonal entries in D (or, respeetively, in U) are positive. 

To do the eomputation with interval arithmetie, we used the mpmath . mpi [19] interval arith¬ 
metie library and the modified sympy [36] funetion SparseMatrix . LDLdecomposit ion 
to eheek positivity of the entries in i9. To do the LU deeomposition with exaet rational arithmetie, 
we used sparse matriees from sage [35] and its built-in LU method. The seeond approaeh gives 
exaet values for the diagonal entries of U, and again, all entries are positive. □ 


8.5. Strategy for computing discrete eigenvalues. The noneonstruetive method deseribed the 
previous seetion already provides two ways to prove Lemma 6.9. Nevertheless, in this seetion, 
we also present a eonstruetive approaeh. As diseussed above, the eigenvalues of the diserete sys¬ 
tem (28) are not available in elosed form. An iterative method will be required to eompute an 
approximation to the eigenvalues. We then eonvert the resulting approximation to bounds for the 
eigenvalues using interval arithmetie. 
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Proposition 8.3. The smallest two mixed eigenvalues for the domains D^f), Duf) lie within the 
ranges shown in Table 2, 3. 

Using interval arithmetic, we first convert to tridiagonal form (with intervals for each of the 
non-zero matrix entries), and then use a bisection method to obtain an approximation to ^Ai/j. 
This method uses the Sturm sequence property [37], and allows us to locate other eigenvalues as 
well. In particular, we also obtain validated approximations to X 2 ,h- 

Then, using interval arithmetic, the Lanczos method [22, 15] on P/^, and a validated nonlinear 
solver we obtain an approximation with error bound for Xi^h- 

For the latter two calculations, we use an interval arithmetic calculation in I NT LAB [31]. Specif¬ 
ically, we begin with the rational matrices K/^, all subsequent matrix operations are done within 
interval arithmetic. 

The use of the bisection method and the Sturm sequence property ensures that we have no miss¬ 
ing eigenvalues in the initial part of the spectrum of (28), and that the first two discrete eigenvalues 
are seperated: 0 < li^h < h,h < h,h- As it standard with this method, the convergence rate is 
slow. In each of the subdomains under consideration, we used this method to verify that (i) the first 
eigenvalue is simple and (ii) the (normalized) spectral gap X 2 ,h — Xi^h is large. This method also 
provides an initial estimate for Ai/j. 

With this initial estimate for Ai^^, we use the Lanczos iteration on the generalized eigenvalue 
problem (26), 

K/iiii/j = Ai^/iM/iUi/i, 

to compute a higher-accuracy approximation to Xi^h and a good approximation to the correspond¬ 
ing eigenvector. Thanks to the (large) spectral gap between the first and the second discrete eigen¬ 
value, we know the first eigenvector will be well-approximated. The specific implementation we 
use is the eigs subroutine in Matlab, which gives as output the eigenpair (Ai/j, ui/i). The val¬ 
idated numerics INTLAB algorithm VerifyEig is then applied to the eigenpair hi./i)- 

The output, using validated interval numerics, is the approximate (interval) eigenpair Xi^h, uj^h- 
The assertion follows from the Banach contraction mapping principle, [32]. Therefore, the true 
eigenvalue Xi^h C {Xi^h ± e} where e is the radius of the eigenvalue inclusion interval. 

In Table 2 and Table 3, we record the midpoints Ai /i, X 2 ^h of the first and second eigenvalues for 
each of the domains. We report the quality of the first eigenvector. 5 is the largest (component¬ 
wise) relative error (radius of interval scaled by midpoint). In Table 4 and Table 5 we record first 
eigenvalue (midpoint of interval) Ai/i, the radius e of this interval, and the quantity Ai— A — L*. 
It is clear that 

Ai,/i - A - L* > e > 0. 

Since the true discrete eigenvalue is provably guaranteed to lie G {Xi^h ± e}, we conclude that 
Xi^h > L*. Consequently, A > 7. 
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Ai,h 
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Table 2. Ai^h = 128^Ai,/i, A 2 ,/i = 128^A2,/i, 5 for the upper domains 


Domain 

Ai,h 

A 2 ,h 

128^5 

0 
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1 

12.27959016 

56.23540471 
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Table 3. Ai 

= 128%,;,, A: 

2,h = 128^A2,/i, 

e for the lower domains 


Domain 

Ai,h 

1282 e 

Ai,;, - 1282e - (12.25 + 3/256) 

0 

12.32808937 

4.14743795e-ll 

6.6370623 le-02 
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12.28049706 

4.44586590e-ll 

1.87783061e-02 
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4.25988134e-ll 

8.96106588e-03 

3 

12.28475539 

4.13429291e-ll 
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4 

12.32978563 

4.97983876e-ll 
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5 

12.27032723 

5.32356381e-ll 

8.60847819e-03 

6 

12.26881872 

3.81223941e-ll 

7.09997069e-03 

7 

12.26612593 

3.64668296e-ll 

4.40717954e-03 

8 

12.28907295 

3.70494746e-ll 

2.73542003e-02 

9 

12.28227434 

4.88817875e-ll 

2.05555934e-02 


Table 4. Ai/i, 128^e, Ai ;il28^e — (12.25 + 3/256), for the upper domains. 
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Domain 
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1 
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2 
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4.30731006e-11 
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3 

12.27157804 

4.69331241e-ll 
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4 

12.32234011 

4.108358106-11 

6.062135766-02 

5 

12.30303087 

4.151168306-11 

4.131211786-02 

6 

12.29807616 

5.089795256-11 

3.635741156-02 

7 

12.30521249 

4.226841106-11 

4.349373886-02 

8 

12.29850781 

4.787459326-11 

3.678906486-02 

9 

12.29425300 

4.569500336-11 

3.253425456-02 


Table 5. 128^e, Ai^/il28^e — (12.25 + 3/256), for the lower domains. 
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